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Abstract — In this article we study the frictionless cooling of 
atoms trapped in a harmonic potential, while minimizing the 
transient energy of the system. We show that in the case of 
unbounded control, this goal is achieved by a singular control, 
which is also the time-minimal solution for a "dual" problem, 
where the energy is held fixed. In addition, we examine briefly 
how the solution is modified when there are bounds on the 
control. The results presented here have a broad range of 
applications, from the cooling of a Bose-Einstein condensate 
confined in a harmonic trap to adiabatic quantum computing 
and finite time thermodynamic processes. 



I. INTRODUCTION 

Frictionless atom cooling in a harmonic trapping potential 
is defined as the problem of changing the harmonic frequency 
of the trap to some lower final value, while keeping the 
populations of the initial and final levels invariant, thus with- 
out generating friction and heating. Achieving this goal in 
minimum time has many important potential applications. For 
example, it can be used to reach extremely low temperatures 
inaccessible by standard cooling techniques [1], to reduce the 
velocity dispersion and collisional shifts for spectroscopy and 
atomic clocks [2|, and in adiabatic quantum computation [3j. 
It is also closely related to the problem of moving in minimum 
time a system between two thermal states, as for example in 
the transition from graphite to diamond [4 |. 

The quest for optimal controls that minimize the neces- 
sary time for the cooling process has an interesting history. 
It was initially proved that minimum transfer time can be 
achieved with bang -bang real frequency controls [4 |. Later, 
it was shown that when the restriction for real frequencies is 
relaxed, allowing the trap to become an expulsive parabolic 
potential at some time intervals, shorter transfer times can be 
obtained \5]. Based on these previous works, we formulated 
frictionless atom cooling as a minimum-time optimal control 
problem, permitting the frequency to take real and imaginary 
values in specified ranges |6|. We showed that the optimal 
solution has again a bang-bang form and used this fact to 
obtain estimates of the minimum transfer times for various 
numbers of switchings |6|. We subsequently fully solved the 
corresponding time-optimal control problem and obtained the 
optimal synthesis |7|. 

Recently, it was pointed out that the energy "cost" of the 
cooling process, more precisely the transient excitation energy, 
can impose limits to the possible speed-up ID. For example, 
for a trap which is harmonic near the ground state, but not 
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for higher energies, large transient energies imply perturbing 
effects of anharmonicities and thus undesired excitations of 
the final state. The problem of minimizing the transient 
energy for a fixed transfer time and with unlimited controls 
was considered in |8|, and an interesting conjugate relation 
between these two quantities was revealed. 

In the present article we examine this problem from a 
control-theoretic point of view. We show that the transient 
energy is minimized by a singular control. We also show that 
this singular control is the time-minimal solution for a "dual" 
problem, where the energy is held fixed, elucidating thus 
the conjugate relation between transient energy and minimum 
time. Finally, we examine briefly how the solution is modified 
when there are bounds on the control. The results presented 
here come to meet several examples of singular solutions for 
optimal control problems on spin systems existing in Quantum 
Control Uterature l|9l, lITOl, ||lT|, lITII. 

II. FORMULATION OF THE PROBLEM IN TERMS 
OF OPTIMAL CONTROL 

The evolution of the wavefunction y/{t,x) of a particle 
in a one-dimensional parabolic trapping potential with time- 
varying frequency co{t) is given by the Schrodinger equation 

^2 fnco^(^t) 



2m dx 



(1) 



where H{t) is the Hamiltonian operator of the system, m is 
the particle mass and h is Planck's constant. Consider the 
time evolution with initial frequency a)(0) = (Bo at f = 
and final frequency (o{T) = o>t < 03q at the final time T. 
This corresponds to a temperature reduction (cooling) by a 
factor cOt/cOq- Frictionless cooling corresponds to a path (X){t) 
between these two values so that the populations of all the 
oscillator levels n — 0,1,2, . . . at t = T are equal to the ones 
at f = 0. More explicitly, if 

CO oo 

V/(0,x) = £ c„(0)^„«o(^)^ y,^T,x) = £ c„{T)^^{x), (2) 

n=0 n=0 

where (x) , ^"^^ (x) are the eigenfunctions of the operators 
H{Q),H{T), respectively, then frictionless coohng corresponds 
to 



|c„(0)|2 = |c„(r)|2, „ = 0,1,2, 



(3) 



It was shown in H, using the Lewis-Riesenfeld invariant 
[fT3l, that frictionless cooling is achieved when co{t) satisfies 
the following Ermakov equation lfT4l 



(4) 



with boundary conditions 

^7(0) = !, ^(0)=0, b{0)^0, 
b{T) = Y, b{T)=0, b{T)=0, 



(5) 
(6) 



and 7= ^coq/cot > 1. Here b{t) is a scaling dimensionless 
function describing the expansion of the wavefunction during 
the cooHng process. When the above conditions are satisfied, 
the nth eigenstate of the initial oscillator at f = evolves 
following the "expanding mode" 



-I [n- 



(2«n!Z7)i/2 



exp 
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where //„ is the Hermite polynomial of degree n, and becomes 
eventually the nth eigenstate of the final trap at f = T, up to 
a global phase factor (independent of the spatial coordinate). 
The instantaneous average energy E„{t) — (^„|//(f)|^„) for 
the nth expanding mode is 



En{t) 



{2n + \)Ti 



co^{t)b^ 



1,2 



and the time-averaged energy is 



En = ]- [ E„{t)dt. 

T Jo 



(8) 



(9) 



If we substitute (O in (|9l) and use the boundary conditions for 
b, we obtain 



(2n 



IcoqT 



i)n 

Jo 



0^ 

1,2 



dt. 



(10) 



We would hke to find (o{t) satisfying the Ermakov equation 
(|4]l and the corresponding boundary conditions (|5} and (|6), 
i.e., the sufficient conditions for frictionless cooling (|3]), while 
minimizing the time-averaged energy (fTHt . for fixed final time 
T . If we set 

b . . ©2(f) 
— , Mf = 



X2 



(11) 



and rescale time according to fnew = ^'^^old' obtain the 
following system of first order differential equations 



on the state variables (xi,X2) are equivalent to those for b 
and b, while the boundary conditions on the control variable 
u lead to the corresponding conditions for b. Note that the 
possibility (0^{t) < (expulsive parabolic potential) for some 
time intervals is permitted [51. Also, for fixed final time T, the 
minimization of cost J corresponds to minimizing the time- 
averaged energy (fl4l i. 

If we omit the boundary conditions on u{t) and solve the 
corresponding problem, the minimum cost that we find is a 
lower bound of the minimum cost for the full problem, where 
these conditions are on. In the following we solve the relaxed 
problem, while the study of the full case will be the subject 
of a subsequent publication. 

Problem 1: Given the system (fTZt . (fT3T l with initial condi- 
tion (.^:i(0),Jt:2(0)) = (1,0) and final condition {xi{T),X2iT)) = 
(7,0), 7 > 1, find the control u{t),0 <t<T, T fixed, that 
minimizes the cost J given in JTSl l. 



III. SINGULAR SOLUTION 



The control Hamiltonian is given by 



1 / 1 

'2 V - X 



He = { xl + ] + AiX2 + A2 ( -Mxi + ■ 



(16) 



where, according to Maximum Principle ifTSl . Aq < is a 
constant and the adjoint variables satisfy the equations 



M — — — +M\u + — 

OX\ A'l \ X^ 

M — — — = — A0X2 — Ai . 
0x2 



(17) 
(18) 



Observe that He is linear in u. The switching function is 
<I> = — A2 (note that xi > due to the repulsive force 1 /x\ at 
x\ = 0). Singular arcs are encountered when A2 = for some 
finite time interval. Note that if also Aq = then from ( fTSl l it is 
Ai = 0, which is forbidden since Maximum Principle requires 
(A(), Ai,A2) ^ IfTSl . So Ao 7^ on a singular ai'c and we set 
Ao = -1. For A2 = 0, (UHi reduce to Ai = -\/x\,X\ =X2, 
respectively. Setting Ai = X2 in ( fTSl l we obtain 



xi 

X2 



-'X2, 



xi 



equivalent to the Ermakov equation Q, and 

(2n + l)hao 



where 



J 



x\- 



dt. 



(12) 
(13) 

(14) 
(15) 



The corresponding optimal control problem is: Given the 
system ([T2|, ([T3]i, with initial condition (jci (0),X2(0)) = (1,0) 
and final condition {xi{T) ^X2{T)) = (7,0), find the control 
m(0,0 <t<T, T fixed, with m(0) = l,u{T) = l/f, that 
minimizes the cost J given in (fTsT l. The boundary conditions 



x\ 

x\ 



(19) 



Since the system (fT2] |. (fTSI l is autonomous, the control Hamil- 
tonian is constant. The equation 



2 1 _ 



(20) 



represents a one-parameter family of singular arcs in state 
space. Combining the adjoint equations we find X2 ~ — 1/xj, 
and using this in ( fT3] l we obtain the expression for the control 
on the singular arc 



= — 



(21) 



A. Unbounded Control 

When the control is unbounded, the state can be shifted 
along the lines of constant x\ instantaneously by the use of 
Dirac impulses in u. Such controls have no effect on the 
performance index J directly, since u{t) does not enter the 
cost function. A typical extremal solution involves an initial 
impulse to move the state to the singular arc (at t = 0+), then 
motion along the singular arc until the line xi = 7 is reached, 
then another impulse to move the state to the target point 
(7,0). The fact that the state must arrive at the target point at 
t = T determines the constant c in ( |20] |. which picks out the 
particular singular arc in the one-parameter family of possible 
arcs. Using ( fT2] i and ( |20] |. one can find the time evolution of 
x\ along the singular arc 



2B 



1, 



(22) 



where B = yT^ + T^— 1, in agreement with the result ob- 
tained in [8 1 using the Euler-Lagrange equation. The constant 
defining the singular arc in (l20t is c= (B/r)2-l. The initial 
impulse moving the state to the singular arc is negative. Now 



X2(r-)=ii(r) = 



f - \Jf + 



(23) 



so X2{T- ) > for r < 7^72 - 1 and X2[T-) < for T > 
7\/ 7^ — 1 ■ The final impulse that drives the state to X2(r+) = 
is positive in the first case and negative in the second. 
Characteristic trajectories in state space for these two cases 
are plotted in Figs. [T] and |2] respectively. It can be easily 
shown t hat dcj dT < for T < j^Jf- - 1 and dcjdT > for 
T > 7^ — 1 . The corresponding minimum value is Cjj^jjj = 
-1/7^ and c e [-1/7^,°°), with 00 for T and c 
for T ^ 00, Note that as the final time T varies, the minimum 
cost J is obtained when //^ = ^ c = ^ T = (7^ - l)/2. 
This case is equivalent to finding the minimum J with the 
final time T free. Returning to the fixed time case, we are 
interested in the time-averaged energy E =2J /T, in units of 
(«+ \/2)h0X), which is given by the following expression 



1 

-f2 



7 



(24) 

and is plotted versus T in Fig. |3] It is shown in Appendix lAl 
that dE jdT < 0, so the time-averaged energy is a decreasing 
function of the final time T. 

There is an alternative way to obtain the above results and 
the optimality of the singular solution. The idea is to define 
a "dual" time-optimal problem and use the results regarding 
singular trajectories of such problems [16|. With this in mind, 
we augment the system (fT2l i. ( fTST l with a state variable X3 
satisfying the differential equation 

1 



X3 



x\ 



(25) 



and the boundary conditions X3(0) —xt,{T) —Q, where E >Q 
is fixed. We now seek a control that drives the state of the 
augmented system from x(0) = (1,0,0) to x(0) = (7,0,0),7> 
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Fig. 1. Singulai' trajectories for various values of the final time T when 
7=10 and the control is unbounded. The trajectories characterized by the 
constant values of the control Hamiltonian c = 0, t = — 1/7^, correspond to 
T = (7^ — l)/2, T = Y\/ 7^ ^ li respectively. The initial impulse, driving the 
starting state (1,0) to the singular arc is negative, while the final impulse, 
driving the state from the singular arc to the target point (7,0), is positive. 
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Fig. 2. Singular trajectory for T > 7^ — 1 . Note that final impulse 
negative. 



1, in minimum time T . Observe that the extra state is defined 
such that the equality 



x\- 



1 , 

— \dt 



(26) 



holds, in other words the time-averaged energy is now fixed 
while the final time is free. The corresponding affine control 
system is i = fix) + Mg(x), where the vector fields are 



L 2 

3 ,^1- 



1 



-E)\ g = (0,-xi,0)' 



The next step is to calculate the following determinants 

Z) = det(^,[/,g],[^,[/,g]]) = -2xt 
D'=det(g,[/,g],[/,[/,g]])=4, 
D"=det(g, [/,§],/) 

where the Lie-bracket is \f ,g\ = {dg/ dx)f — {d f/dx)g. The 
singular control for the augmented time-optimal problem is 
given by = D' /D = 2/x\ |[16|, and is the same as in ^2\\ . 
As a consequence, the projection of the singular arc on the 
xiX2-plane is given by (|20| | and T is determined implicitly 
by (l24l i. since E is fixed. We now show that this solution 
corresponds to a minimum. We first show that the singular 



80 




Fig. 3. Time-averaged energy E, in units of («+ \/2)%(Oo, plotted as a Fig. 4. Trajectories con'esponding to tlie (bounded) control policy \i5\ . for 
function of the final time T, for 7=10 and unbounded control. Note that various values of the final time T > Ty^^^ and 7 = 10. 
dE/dT < 0. 



arc is hyperbolic. The criterion for this is llT6l 

DD" >O^E+xl-\=E + c>0. 



(27) 



Using ( |24] | and the expression for c given above we find 
E + c = 2f{T)/T^, where the function f{T) > IS given in 
Appendix |A] So the singular arc is hyperbolic. In this case, 
we know that it corresponds to the minimum time solution, 
up to the first conjugate time ti^ [16|. In order to find fi^ 
we introduce the vector field S = f + u^g, that generates the 
singular trajectory, and the Jacobi field V, which is the solution 
of the variational equation 

5x — -^5x, 
ox 

with the initial condition 5x{0) ^ g{0) = (0,-1,0)^. The first 
conjugate time is the first time such that the fields V,g are 
coUinear, i.e., V(fk) || (0, 1,0)^. The first two equations of the 
variational system are 

5xi = 5x2, 
5x2 = -j5xi, 

A 

with the initial conditions 5xi{Q) ~Q,5x2{Q) = —1. It is not 
hard to see that 5xi{t) < for f > 0, so V cannot be aligned 
with g and there is no conjugate point. As a result, the singular 
solution is time-optimal for the augmented system with fixed 
time-averaged energy. This elucidates the conjugate relation 
between energy and time. 

B. Bounded Control 

In practice there are limits on the control amplitude that set 
a tighter lower bound on the transfer time. To fix ideas, we 
consider the case \u\ < 1 for the two-dimensional system (fT2l l. 
( fT3] l. We have shown in our previous work [7] that in this 
case the minimum necessary time to transfer the initial state 
(1,0) to the target state (7,0), 7 > 1, is obtained following the 
bang-bang strategy 



t < Ti 

Ti<t<Ti+T2 



(28) 



From (fT2] i. (fTSl l we find the trajectory for t <Ti 

xl~x] + ^^0, 
x\ 

and the trajectory for Ti <t <Ti +T2 

4 



The xi coordinate of their common point at t = Ti is 



xi(ri) = 



if + i 



If 



(29) 



(30) 



(31) 



Integrating (fT2l i. ( fTsl l from the initial point up to the common 
point, for duration t ~T\ and m = — 1 , we find 



Xl 



(ri) = v/cosh(2ri). 



(32) 



Integrating (fT2] |. dTsl l from the final point back to the common 
point, for duration t = T2 and m = 1 , we find 



^1 (ri) = J ^ [f + ^)+\{f-Y2] '^°^(2r2) (33) 



From dSTT l. ( |32] | and ( l33T l we obtain 



1 



■ cosh 



272 



K 

4 



(34) 



We see that in the presence of control bounds, problem [T] is 
meaningful for final times T > T^jjj- For simplicity, we will 
restrict the rest of the discussion to the bang-singular-bang 
case 



(35) 



-1, t<Tl 
u{t)^{ II x\, Tl<t<Tl + Tl 

1, T! + T;<t<T! + T; + n 



Under this control policy, the initial and final arcs (129) and 
( |30l l are joined by a singular arc (|20t . as it is depicted in Fig. 
|4] We expect that at least for the most interesting case, where 
T is close to the pulse sequence ( [35] ) is the optimal one. 

Note that the joint point for T — T^jfj lies above the singular 
arc X2 = l /xi corresponding to c = 0, as shown in Fig. |4] when 
7**-67^+l>0^7> \/3+2V2 w 1.554. Throughout the 
text we use the value 7= 10. 
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Fig. 5. Time-averaged energy corresponding to the (bounded) control policy 
)35K as a function of the final time T > T^^^, for 7= 10 (dashed line). The 
con'esponding energy for the case of unbounded control and the same target 
point is also plotted (solid line). 



Using (|29T l. (|20] | we find for the first junction point 



(36) 



while from (|30] |, (|20] | we find for the second junction point 



-cY^ + f+l + ^y{cY^-f-l)^-Sf 
2f ■ 



(37) 

If we limit our analysis to the cases where the constant deter- 
mining the singular arc satisfies c > 0, which also correspond 
to shorter transfer times, then 



< 



< 1, 



(38) 



so the control remains within the allowed bounds along the 
singular arc. By integrating the equations of motion we obtain 
the time spent on each arc 



T(ic)^-cosh-'xliTl), 



1 + cx\{T[ + T{) + a/1 + cx\{T[) 



(39) 
(40) 



The constant c which determines the singular arc can be found 
from the transcendental equation 



T[{c) + T{{c) + Ti{c) = T. 



(42) 



Having determined c, we can calculate the time-averaged en- 
ergy corresponding to the control sequence (l35i for a specific 
final time. The result is plotted in Fig. |5] along with the time- 
averaged energy that derived previously for the unbounded 
control case. As expected, when the control is bounded, a 
larger amount of energy is necessary to achieve the same 
transfer for times close to T^jjj- Equivalently, for a fixed level 
of time-averaged energy, a longer time is needed to reach the 
same target point. 



IV. CONCLUSION AND FUTURE WORK 

In this paper we examined from a control-theoretic view- 
point the frictionless cooling of atoms trapped in a harmonic 
potential, while minimizing the transient energy of the system, 
a problem that was first considered in [S]. We showed that 
the transient energy is minimized by a singular control. We 
also showed that this singular control is the time-minimal 
solution for a "dual" problem, where the energy is held fixed, 
highlighting the conjugate relation between transient energy 
and minimum time. In addition, we examined briefly how the 
solution is modified when there are bounds on the control. 

Possible future work could include the detailed examination 
of the bounded control case, that was only slightly touched 
here, as well as the incorporation in the analysis of additional 
restrictions on the controls reflecting experimental limitations 
ifTTl . The complexity of the resulting optimal control prob- 
lems, which may increase the difficulty of the analytical study, 
can be overcome by using a powerful state of the art numerical 
optimization method based on pseudospectral approximations 
PJSI, [19|. 

The results presented here can be immediately extended to 
the frictionless cooling of a two-dimensional Bose-Einstein 
condensate confined in a parabolic trapping potential |20|, 
1,21]. The above techniques are not restricted to atom cooling 
but are applicable to areas as diverse as adiabatic quantum 
computing [[3] and finite time thermodynamic processes [|4j. 
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Appendix A 

We show that dE/dT < 0. From (EUl we obtain dE/dT = 
-2/(^)/^^ where 



f{T) = Y^ + l-2^fTf^+Tln ^V/+ ^43) 

7 



It is /(O) = (7- 1)2 > and df/dT = g{T), where 



It is g{0)=0 and 



dT 



T^-{f + T^)-y^>0. 



(44) 



(45) 



So g{T) > ^ f{T) >0^ dE„/dT < 0. 
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